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g : ABSTRACT 

o 



We describe observations of the nearby (cz = 11,487 km s~^) cluster of 
galaxies Abell 576 beyond the virial radius and into the infall region where 
^ I galaxies are on their first or second pass through the cluster. Using 1057 



redshifts, we use the infall pattern in redshift space to determine the mass profile 
of A576 to a radius of ~ 4 h~^Mpc. This mass estimation technique makes no 
assumptions about the equilibrium state of the cluster. Within ~ 1 h^^Mpc, 



H ' the mass profile we derive exceeds that determined from X-ray observations by 

a factor of ~ 2.5. At ~ 2.5 h Mpc, however, the mass profile agrees with virial 
mass estimates. Our mass profile is consistent with a Navarro, Frenk, & White 
(1997) or Hernquist (1990) profile, but it is inconsistent with an isothermal 
sphere. 

R-band images of a 3° x 3° region centered on the cluster allow an 
independent determination of the cluster light profile. We calculate the 
integrated mass-to-light ratio as a function of cluster radius; it decreases 
smoothly from the core to M/Lr ~ 300/i at ~ 4 /i-^Mpc. The differential 
M/Lji profile decreases more steeply; we find 5M/5Lr ~ lOO/i at ~ 4 /i~^Mpc, 
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in good agreement with the mass-to-hght ratios of individual galaxies. If the 
behavior of M/Lr in A576 is general, VL^ <^ 0.4 at 95% confidence. 

For a Hernquist model, the best-fit mass profiles differ from the observed 
surface number density of galaxies; the galaxies have a larger scale radius than 
the mass. This result is consistent with the centrally peaked M/Lr profile. 
Similarly, the scale radius of the light profile is larger than that of the mass 
profile. We discuss some potential systematic effects; none can easily reconcile 
our results with a constant mass-to-light ratio. 



Subject headings: galaxies: clusters: individual (A576) — dark matter — 
galaxies: kinematics and dynamics — galaxies: photometry — cosmology: 
observations 



Introduction 



Clusters of galaxies are important probes of the distribution of both matter and 
light on intermediate scales (0.1 — 10 /i~^Mpc). They are also interesting laboratories 
for studying the effects of local environment on cluster galaxies. Most studies of clusters 
concentrate on the central regions where the cluster is probably in equilibrium. Studies 
of the galaxy distribution on larger scales tend to focus either on general properties of 
large-scale structure (e.g., [deLapparent, Geller, fc Huchra 1986| ; pressler et al. 19871) or 
on individual superclusters (e.g., Giovanelli & Haynes 1985). There are relatively few 
examinations ( [Regos fc Geller 1989 ; Lilje fc Lahav 1991 ; [van Haarlem et al. 1993| ; 



van 



Haarlem fc van de Weygaert 1993t Praton fc Schneider 199^ ; piaferio fc Geller 1997| ; [Vedel 
fc Hartwick 1998| ; pilingson et al. 1999|) of the infall regions of clusters. In these regions, the 
galaxies are falling into the gravitational potential well of the cluster, but they have not yet 
reached equilibrium. Many, perhaps most, of the galaxies in this region are on their first 
orbit of the cluster. They populate a regime between that of relaxed cluster cores and the 
surrounding large-scale structure where the transition from linear to non-linear clustering 



occurs. 



Here we use 1057 redshifts and photometry of 2118 galaxies in the infall region of Abell 
576 {cz = 11,487 km s~^) to address an unresolved problem in astrophysics: the relative 
distribution of mass and light on large scales. Zwicky (1933; 1937) originally found that the 
mass of the Coma cluster greatly exceeds the sum of the masses of the stars. More recently, 
calculations of mass-to-light ratios for galaxy clusters yield values of several hundred in solar 
units ( pressler 1975| ; [b'aber & Gallagher 1979| ; |Adami et al. 1998 b| ). Girardi et al. (2000) 
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calculated mass-to-light ratios for a large sample of nearby clusters within the virial radius 
and obtained typical values of M/ Lb^ ~ 220 — 250/i (all M/L values are in solar units, i.e., 
Mq/Lq). David, Jones, & Forman (1995) find M/Ly ~ 200 — 300/i for seven groups and 
clusters using masses calculated from the observed X-ray emission, and the CNOC survey 



finds M/ Lr ~ 290 ± 60/;, for a sample of distant clusters QCarlberg et al. 1996|) 



The mass density parameter VLm of the universe can be estimated by assuming that 
the universal mass-to-light ratio is equal to the ratio in rich clusters of galaxies (e.g., 
Carlberg et al. 1996). There have been few attempts to determine mass-to-light ratios 
directly at larger radii ( |Mohr fc Wegner 1997| ; [Small et al. 199^ ; [Kaiser et al. 20001) , where 



clusters are in neither hydrostatic nor virial equilibrium. Neither X-ray observations nor 
virial analysis provide accurate mass determinations at these large radii. Two methods 
with particular promise are weak gravitational lensing ( [Kaiser et al. 2000] ) and kinematics 



of the infall region (Diaferio & Geller 1997, hereafter DG; Diaferio 1999, hereafter D99). 
Kaiser et al. analyzed the weak lensing signal from a supercluster at z ~ 0.4 and found no 
significant evidence for variations in mass-to-light ratios on scales less than ~ 6 /i^^Mpc. 
Metzler et al. (1999), however, show that lensing by intervening filamentary structures 
probably associated with clusters can result in significant overestimates of cluster masses. 
Geller, Diaferio, & Kurtz (1999, hereafter GDK) applied the kinematic method of DG to 
the infall region of the Coma cluster. They successfully reproduced the X-ray derived mass 
profile and extended direct determinations of the mass profile to ~ 10 /i^^Mpc. Adequate 
photometric data necessary to compute the mass-to-light profile for this large region around 
Coma are not yet available. 

Here we apply the method of DG to A576. In redshift space, the infall regions 
of clusters form a characteristic trumpet-shaped pattern. These caustics arise because 
galaxies fall into the cluster as the cluster potential overwhelms the Hubble flow ( [Kaiser[ 
1987[ ; [Regos fc Geller 1989[ ). Under simple spherical infall, the galaxy phase space density 



becomes infinite at the caustics. DG analyzed the dynamics of infall regions with numerical 
simulations and found that in the outskirts of clusters, random motions due to substructure 
and non-radial motions make a substantial contribution to the amplitude of the caustics 
which delineate the infall regions. DG showed that the amplitude of the caustics is a 
measure of the escape velocity from the cluster; identification of the caustics therefore 
allows a determination of the mass profile of the cluster on scales ^ 10/i~^Mpc. 

DG show that nonparametric measurements of caustics would yield cluster mass 
profiles accurate to ~30% on scales 1 — 10 /i~^Mpc, if (1) the redshift space coordinates of 
the dark matter particles were measurable, and (2) the cluster mass within the virial radius 
were known exactly. More realistically, by using simulated catalogs of galaxies formed 
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and evolved using semi-analytic procedures within the dark matter halos of dissipationless 
N-body simulations ( |Kautl:'mann et al. 1999a| ), D99 shows that the identification of caustics 



in the realistic redshift diagrams of clusters recovers their mass profiles within a factor 
of 2 to several Megaparsecs from the cluster center without a separate determination 
of the central mass. When combined with wide-field photometry, this approach allows 
a determination of the mass-to-light ratio on large scales which is independent of the 
assumption that light traces mass. This method assumes only that galaxies trace the 
velocity field. Indeed, simulations suggest that velocity bias, if any, is very weak on both 
linear and non-linear scales ([Kaufi'mann et al. 19994 Piaferio et al. 1999D . Vedel & Hartwick 



(1998) used simulations to explore an alternative parametric maximum likelihood analysis 
of the infall region. Their technique requires assumptions about the functional forms of the 
density profile and the velocity dispersion profile. 

Here, we analyze the caustics within a ~ 4 h~^Mpc radius to determine the mass profile 
of Abell 576, an Abell richness class 1 cluster ([Abell 1958| ) at a redshift of z = 0.0383. Mohr 
et al. (1996, hereafter M96) extensively studied the inner square degree of this cluster. 
Using photometric observations of a 3° x 3° region centered on the cluster, we determine 
the mass-to-light profile within this range. The integrated mass-to-light ratio smoothly 
decreases with radius out to ~ 4 h~^Mpc. Remarkably, the differential mass-to-light ratio 
decreases steeply to values typical of individual galaxy halos. 

We describe the photometric and spectroscopic observations in § 2. In § 3, we 
determine the amplitude of the caustics and calculate the mass profile. We analyze the 
X-ray observations in § 4 and compare the X-ray mass profile to the infall mass profile. We 
discuss the contribution of background galaxies in § 5. We determine the surface number 
density profile in § 6 and compare it to the mass distribution. We test for luminosity 
segregation in § 7 and calculate the light profile in § 8. We compare the mass and light 
profiles of the cluster in §9. We discuss possible systematic effects in §10 and conclude in 
^11. 



2. Data 
2.1. Images 

We obtained ten 200-second R-band images of A576 with the MOSAIC camera on the 
0.9-m telescope at Kitt Peak on 1999 February 13-14; both nights were photometric. The 
MOSAIC camera consists of 8 CCD chips, each with a field of view of 15' x 30'. Our mosaic 
thus covers approximately a 3° x 3° region centered on A576. Adjacent images overlap by 
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I'O. We imaged the central square degree on both nights; the 14 February image is offset 
I'O N and I'.O E from the 13 February image. We reduced the images using standard IRAF 
procedures in the mscred package. We use SExtractor ( [Bertin fc Arnouts 1996D to locate 



sources in the images and calculate magnitudes. SExtractor divides the image into segments 
and assigns pixels to individual objects to deal with crowded images. Analysis of simulated 
CCD images show that SExtractor accurately recovers total magnitudes of faint galaxies 
( Pertin fc Arnouts 1996|) . We use the MAG_BEST magnitudes, which are equivalent to 
MAG_AUTO, an aperture magnitude, unless an object has one or more close neighbors 
likely to contaminate the flux by more than 10%, in which case MAG_BEST=MAG_ISOC, 
an isophotal magnitude with a correction based on a Gaussian light profile. We classify 
all objects with CLASS_STAR< 0.6 as galaxy candidates; we then visually inspect these 
to eliminate binary stars and artifacts. We calibrate the photometry with observations of 
M67 ( [Chevalier fc llovaisky 1991| ; |Anupama et al. 1994|) and Landolt (1992) fields. We 



obtain the color correction term from the data set for the entire observing run and the 
extinction coefficient from fits to data on all chips for each night. To allow for possible 
chip-to-chip sensitivity variations, we calculate the zero point separately for each chip 
( IBrown et al. 200q) . 

We use the IRAF package apphot to test the accuracy of the SExtractor magnitudes. 
Using curves of growth on the images of a few dozen galaxies covering the range of 
magnitudes rriji = 14 to rriji = 18, we estimate their total magnitudes within large apertures 
(50-150"). We find excellent agreement (^ 0.1 mag) between the IRAF total aperture 
magnitudes and the SExtractor MAG_BEST magnitudes. 

We restrict our analyses to objects with MAG_BEST<18.0, approximately the limit of 
the classification. We use the offset images of the central region to estimate the consistency 
of our magnitudes; the distribution of m2 — mi (mi and m2 are the magnitudes calculated 
from the first and second night respectively) for galaxies with m^j <18.0 is well represented 
by a Gaussian with a zero point difference of -0.002 mag and (J(m2-mi) ~ 0.09 mag. The 
scatter in the magnitude determination increases significantly with apparent magnitude 
(Figure p. We thus estimate that our magnitudes are internally consistent to about 0.06 
mag (cr^^ ~ (7^2 ~ 0"(m2-mi)^/2). In Table 1, we list the R band magnitudes and their 
uncertainties (quadrature sums of the internal SExtractor errors and the 0.06 mag scatter 
from galaxies with multiple observations) for galaxies with redshifts from FLWO (§2.2). We 
also include estimates of Galactic extinction in the R band {An) based on the dust maps of 
Schlegel, Finkbeiner, & Davis (1998). Mohr et al. (2000) will list redshifts collected by JJM 
and GW. Magnitudes for these galaxies as well as those without redshifts are available upon 
request from the authors. In the entire survey region, we identify 825 and 2118 galaxies 
with niji <17.0 and m^j <18.0 respectively. 
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We compare our photometry with that of M96. Figure ^ shows the difference in 
magnitude between the two studies as a function of our magnitudes. The comparison 
suggests that our magnitudes are systematically brighter by ~ 0.1 mag. Further analysis, 
however, shows that the differences between our magnitudes and M96 can be attributed to 
differences in the reduction packages (SExtractor versus FOCAS). For example, FOCAS 
uses a global estimate of background counts while SExtractor uses a local estimate. A 
complete comparison of the differences between these packages is beyond the scope of this 
work. A reanalysis of the M96 data using SExtractor produces excellent agreement (Figure 
0b) between the two independent photometric data sets ( |Mohr et al. 200(]| ). 

Out of 823 galaxies with redshifts selected without reference to these images (§2.2), 
SExtractor misses 28. Visual inspection reveals that 17 of these galaxies are in the halos 
of bright stars; the rest lie in chip gaps. For the 17 galaxies in stellar halos, we perform 
photometry with apphot. Unfortunately, three of these galaxies are located within ~ 1'.5 
of the X-ray center of the cluster (a saturated star is I'.l from the X-ray center), which 
increases the uncertainty of the central light. We make no correction for the regions of sky 
covered by saturated stars; this approach may lead to a slight underestimate in the light 
profile. 

2.2. Spectroscopy 

We used the FAST spectrograph ([Fabricant et al. 1998|) on the 1.5-m Tillinghast 
telescope of the Fred Lawrence Whipple Observatory (FLWO) to obtain 529 spectra of 
galaxies within 4° (f» 8 h~^Mpc) of the center of A576. FAST is a long-slit spectrograph 
with a CCD detector. Integration times were typically 4-20 minutes, with spectral 
resolution of 6-8 A. We analyzed the spectra at the Telescope Data Center at the 
Smithsonian Astrophysical Observatory using the xcsao and emsao tasks, which are 
standard cross-correlation fitting routines ( [Kurtz fc Mink 1998|) available in the IRAF 
package r2rvsao. In this technique, template galaxy spectra are cross-correlated with the 
observed log-wavelength binned spectra to determine the redshift which provides the best 
fit for a particular template. Cross-correlation fits with R values of ^ 4.0 are acceptable; 
most spectra have larger R values (the goodness of fit increases with R value). 

We observed infall galaxy candidates in two campaigns; these occurred before and after 
obtaining the MOSAIC images described above. To select candidates for the first campaign 
(1999 January- February), we obtained galaxy positions from digital scans of the POSS I 
plates, using a method developed by Daniel Koranyi (1999) to discriminate between stars 
and galaxies. We extracted sources from the plate scans with SExtractor. We included all 
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objects with staricity parameter <0.6 and IS07-to-IS02 isophotal ratio greater than 0.25. 
Comparison with visual inspection suggests that ^ 1% of galaxies are misclassified as stars 
using this criterion. We visually inspected the remaining objects to eliminate stars from the 
sample. For galaxies with mn ^ 15.3, we measured redshifts in magnitude order. Because 
the magnitude order was taken from the plate scans, this ordering is not strictly accurate, 
and we could not obtain redshifts for some low surface brightness galaxies. Out of 300 
candidates, this campaign yielded 293 redshifts (11 candidates were stars, plate flaws, or 
unobservable; 4 galaxies were serendipitously observed). 

In the second campaign (1999 October-2000 February), we selected candidates from 
the MOSAIC images. Prior to the second campaign, our redshift sample was ~ 95% 
complete for niR < 15.3. In the second campaign, we observed 236 galaxies in magnitude 
order as determined from the MOSAIC images. We list the redshifts and uir magnitudes 
of the galaxies from both campaigns in Table 1. Table 2 lists the positions and redshifts of 
galaxies outside the MOSAIC images. 

Two of us (JJM and GW) obtained 528 redshifts of galaxies in the central 2° x 2° of 
this region for a separate Jeans' analysis of the central region of A576 ( |Mohr et al. 2000| ). 



281 of these redshifts are from M96; the remainder will be published in Mohr et al. (2000). 
These redshifts were measured at the Decaspec ([Fabricant fc Hertz 1990| ) at the 2.4-m 
MDM telescope on Kitt Peak and Hydra, the multifiber spectrograph on the WIYN 
telescope on Kitt Peak. Figure ^ shows the fraction of galaxies with redshifts as a function 
of R-band magnitude. We divide the sample into the region mostly covered by JJM and 
GW (projected radius Rjj < 1°) and the region mostly covered by this study {Rp > 1°). 
Within the 3° x 3° region covered by our images, the total sample of 817 galaxies is 100% 
(99.7%) complete to a limiting magnitude of = 16.2(16.5); the latter is the completeness 
limit for the photometric region. 



3. Galaxy Infall Method and the Mass Profile 

We briefly review the method DC and D99 develop to estimate the mass profile of 
a galaxy cluster by identifying its caustics in redshift space. The method assumes that 
clusters form through hierarchical clustering and requires only galaxy redshifts and positions 
on the sky. The amplitude A{r) of the caustics is half of the distance between the upper 
and lower caustics in redshift space. Assuming spherical symmetry, A{r) is related to the 
cluster gravitational potential 0(r) by 
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where /3(r) is the velocity anisotropy parameter. DG show that the mass of a spherical shell 
of radii [ro,r] within the infall region is given by the integral of the square of the amplitude 



where Fp ^ 0.5 is a filling factor with a numerical value estimated from simulations. 
Variations in lead to some systematic uncertainty in the derived mass profile. Our mass 
profile extends to ~ 3r2oo (the average density within rs is 6 times the critical density, 
§4.1); within this radius, varies by <^ 15% in simulations (see D99 for a more detailed 
discussion) . 

Operationally, we identify the caustics as curves which delineate a significant drop 
in the phase space density of galaxies in the projected radius-redshift diagram. Galaxies 
outside the caustics are also outside the turnaround radius. For a spherically symmetric 
system, taking an azimuthal average amplifies the signal of the caustics in redshift space 
and smooths over small-scale substructures. D99 described this method in detail and 
showed that, when applied to simulated clusters with galaxies modelled with semi-analytic 
techniques, it recovers the actual mass profiles within a factor of 2 to 5-10 h~^Mpc from the 
cluster center. D99 describes some potential systematic effects including projection effects 
and variation in the galaxy orbit distribution (3{Rp). 

Simulations are necessary to estimate the uncertainties due to projection effects and 
deviations from spherical symmetry. In the simulations of D99, the degree of definition of 
the caustics depends on the underlying cosmology; in simulations, caustics are better defined 
in a low-density universe than a closed, matter-dominated universe (D99). Surprisingly, the 
caustics of Coma, A576, and several other clusters ( [Rines et al. 2000|) are generally better 
defined than those of the simulated clusters. Thus, the uncertainties estimated from these 
simulations might be overestimates. 

We apply the technique of D99 to our A576 survey to determine the spatial and velocity 
center of the region as well as the location of the caustics. In this technique, we determine 
the center of the system from a hierarchical cluster analysis. The center thus derived, Uopt = 
7:21:31.96, 5opt =+55:45:20.6 (J2000), cz = 11,487 km s'^ lies 50" (28 h-^kpc) from the 
X-ray center (M96, from Einstein IPC data; see also § 4 for ASCA data). Our measurement 
of the central velocity of the infall region agrees well with previous estimates (Struble & 
Rood 1991, M96). 

Figure ^ shows the projected radius from the cluster center versus redshift for galaxies 
within cziim = 4000 km s~^ of the cluster center. This range of redshifts includes all galaxies 
which may be members of the infall region. Solid lines indicate the caustics determined 



Air) 




(2) 
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using the method of D99 based on a multidimensional adaptive kernel method ( [Silverman 



|1986| ; [Pisani 1993| ; [Pisani 1996| ). This technique has been applied to numerical simulations 



(D99) as well as to the Coma cluster (GDK). As discussed in D99, it is necessary to rescale 
hy and hr-, the velocity and radial smoothing lengths so that spherical smoothing windows 
can be used. With an appropriate choice of this scaling relation q = h^/hr, the location 
of the caustics should be insensitive to small changes in q. In our case, g = 25 satisfies 
this criterion. Following D99, we define the caustic amplitude A{r) as the minimum of the 
upper and lower amplitude estimates. This prescription is identical to averaging the two 
estimates for an isolated, spherically symmetric system; taking the minimum reduces the 
sensitivity to massive substructure and contamination. 



We note that in A576, the location of the caustics is sensitive to substructure at 
1.5 — 2.2/;,^^Mpc. The caustic amplitude decreases sharply to ~ 800 km s~^ somewhere 
in this region, though the radius of the decrease is sensitive to the smoothing parameter 
q. Figure |^ shows the caustics determined by setting q = 10, 25, and 50. Small values of 
g ~ 10 seem to oversmooth the caustics and the sharp decrease occurs at ~ 2.2/i~^Mpc, 
whereas the sharp decrease occurs at ~ 1.5/i~^Mpc for larger values of g = 25 — 50. The 
amplitude of the caustics is stable both at radii smaller than ~ 1.5/i~^Mpc and at radii 
larger than ~ 2.2/i-^Mpc. D99 gives the prescription for estimating the uncertainties in 
the caustic amplitude; this prescription reflects the scatter due to projection effects in the 
simulations. We show these la uncertainties in Figure ^. 

We define membership of the infall region from the caustics; hereafter, galaxies outside 
these caustics are interlopers. Of the 497 galaxies within czum of the velocity center of 
A576, 368 are within the infall region. Figure ^ shows the distribution of these galaxies on 
the sky. There is a noticeable deficit of infalling galaxies NW of the cluster center. 



3.1. Comparison of Mass Profile to Models 

We compare our results to Navarro, Frenk, & White (1997, hereafter [NFW| ), Hernquist 
(1990), and singular isothermal sphere mass profiles. The Hernquist profile is an analytic 
form proposed as a model of elliptical galaxies and bulges (Hernquist 1990). The "universal 
density profile" of NFW accurately models the mass profiles of dark matter halos in a 
variety of cosmological simulations (|NFW]). Other simulations suggest that the NFW profile 



may not be accurate at small radii (e.g., [Moore et al. 199^ ), but these differences are 



unimportant on the scales we probe here. These mass profiles are: 

\ — ]i (3) 
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Muerni^) — CRern-, ; TTT, (4) 

(r + a)^ 

Miso{r) = QsoT (5) 

respectively, where a is the characteristic radius and C is a normahzation factor. These 
forms of the mass profiles minimize the correlation between the parameters. For the 
Hernquist profile, the mass Mc within a is CnemO'/'^i and the total mass of the system 
is Mtot = CuernO'- From Equation |^, a singular isothermal sphere mass profile produces 
caustics with constant amplitudes. However, Figure ^ shows that the amplitude of the 
caustics decreases with radius. 

We fit these models to the observed mass profile by minimizing x^; Table |^ lists the 
results. The measures of the cumulative mass profiles are not independent; thus, the values 
of are indicative and are only meaningful when compared with each other. The best-fit 
parameters are insensitive to the inclusion or exclusion of the three data points where the 
caustics are unstable (1.5 < Rp < 2.2/i~^Mpc). Figure |^ shows the three profiles which best 
fit the infall mass profile and Figure |^ shows the mass density profile. The latter display 
has the benefit that the data points are largely indepedendent of one another. The density 
of the cluster varies by five orders of magnitude across our sample; the overall agreement 
between these profiles and the data is remarkable. The isothermal sphere profile is strongly 
excluded. The Hernquist profile yields a better than the NFW profile, but both are 
acceptable. This conclusion agrees with GDK, who found that the mass profile of Coma is 
much better described by an NFW profile than an isothermal sphere. 

From the infall mass profile, we can derive the values of overdensity radii rg directly. 
The mass Ms contained within the overdensity radius rs exceeds the critical density by a 
factor of 6. Two of the most commonly used radii are r2oo and t^qq. The values of these 
radii are usually determined indirectly. For instance, Carlberg et al. (1997, hereafter CYE) 
define r2oo using the velocity dispersion a of a cluster. Evrard, Metzler, & Navarro (1996) 
use simulations to calibrate a correlation between rsoo and the emission-weighted X-ray 
temperature Tx- We estimate rsoo = 0.96 ± 0.05/i~^Mpc and r2oo = 1-42 ± 0.07/i^^Mpc. 
A576 has an emission- weighted temperature of Tx = 3.77 keV (see §5), which yields an 
estimate of rsoo = 0.76 ± 0.16/i~^Mpc using the Evrard et al. estimator, in agreement with 
our result. 



3.2. Velocity Dispersion Profile 



Many recent papers analyze the velocity dispersion profiles of clusters (e.g., [I'adda 
et al. 1996| ; |den Hartog fc Katgert 1996| ; M96). When combined with the galaxy number 
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density profile in the Jeans equation, the velocity dispersion profile can provide an estimate 
of the mass profile; Mohr et al. (2000) will perform this analysis for A576. Figure ^ shows 
the velocity dispersion profile of A576 where we compute the dispersions in bins of 25 
galaxies. We also display the cumulative projected velocity dispersion profile (Tp{< Rp) 
(calculated from all galaxies inside Rp). The profile is centrally peaked and would probably 
be classified as 'Peaked' by den Hartog & Katgert (1996). 

Several authors ( [Fadda et al. 1996t CYE) suggest that the velocity dispersion of a 
cluster is best estimated by the asymptotic value of the cumulative projected velocity 
dispersion. For A576, (Tp{< Rp) decreases monotonically with radius for Rp ^ 0.5h~^Mpc, 
suggesting that an asymptotic value of (Jp{< Rp) may not exist. At the largest radius 
we study, o'p{< Rp) ~ 800 km s^^. This velocity dispersion is significantly smaller than 
ap = 977tlf^ km s"^ (M96) or ap = 9Utll km s"^ (Girardi et al. 1998). 

The CYE estimate of r2oo is sensitive to the definition of ap. We obtain 
^"200 = 1-69 ± 0.17/i~^Mpc using the M96 value of ap, 1.5a larger than our estimate 
of rsoo = 1-42 ± 0.07/i-iMpc (§3.1). Taking the value of ap{< Rp) ^ 800 km s"^ from the 
limit of our survey, r2oo = 1-36 ± 0.15/i^^Mpc, in agreement with our infall estimate. The 
sensitivity of the CYE estimator to the aperture used to measure ap may affect many of 
the results from studies of CNOC clusters because the properties of the composite CNOC 
cluster depend on the estimates of r2oo for individual clusters. 

Based on the infall mass profile, we can predict the velocity dispersion profile with an 
assumption about the distribution of galaxy orbits. We assume that orbits are isotropic 
at all radii (/5 = 0), and calculate the velocity dispersion profiles for the Hernquist and 
NEW models ( [Hernquist 1990| , Equation 10; NFW, Equation 3). These predicted velocity 
dispersion profiles agree with the observed profile (Figure The largest differences are at 
the radii where the caustic amplitude is unstable. Without fitting any parameters, we find 
= 42 and 59 for 14 degrees of freedom for the Hernquist and NFW profiles respectively. 
This result provides a consistency check on the infall mass profile and confirms that a 
Hernquist profile models the data better than an NFW profile. 

3.3. Alternative Kinematic Mass Estimators 

The two most commonly used kinematic mass estimators for clusters are the virial 
mass estimator and the projected mass estimator ( [Heisler, Tremaine, fc Bahcall 19"85|) . 
These estimators both assume that clusters are relaxed systems. Numerical simulations, 
however, suggest that both of these mass estimators typically overestimate the true mass 
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profile of a relaxed cluster by ~ 20% (|Aceves fc Perea 1999|) . 



If clusters are not relaxed, these estimators may provide more substantial overestimates 
of the mass. Figure |^ shows the virial mass estimates of A576 calculated by Girardi et 
al. (1998) from the M96 data. They report both the standard virial mass and a corrected 
virial mass which includes a correction for the surface pressure term estimated from the 
galaxy distribution. This correction assumes that light traces mass, or more precisely, that 
the number density of galaxies traces mass. Girardi et al. also use the galaxy number 
density profile to estimate the mass of A576 at small radius for comparison with X-ray 
estimates. 

M96 show that their data allow a wide range of masses (0.6 — > 1.5 x IO^^H^^Mq) 
depending on the mass estimator, the magnitude cutoff, and the definition of cluster 
members. In particular, emission dominated galaxies have a larger velocity dispersion than 
absorption dominated galaxies. The virial mass of emission line galaxies is a factor of ~ 2 
larger than the virial mass of absorption line galaxies (M96; [Uarlberg at al. 1997|) . Note, 
however, that when we combine the galaxy number density and the velocity dispersion 
profiles in the Jeans equation, the two subsamples should yield consistent mass profiles 
(e.g., [Uarlberg at al. 1997 ). M96 conclude that their data are insufficient to constrain the 



mass of A576 well. The M96 mass range encloses the mass estimates and uncertainties 
given by Girardi et al. (which make no correction for galaxy populations) as well as our 
infall mass estimate. 

Restricting our analysis to galaxies within r2oo = 1.42/i~^Mpc (§3.1), we use the 
projected mass estimator to estimate M = (11.1 ± 2.2) x 10^^ h~^MQ within r2oo assuming 
isotropic orbits; applying the virial theorem yields M = (10.4 ± 2.0) x 10^^ /i~^Mq within 
ryir = 1.15 ± 0.2/i~^Mpc. We display all of these estimates in Figure |^ All of the mass 
estimates at large radius exceed the infall estimate, although the corrected estimate by 
Girardi et al. is consistent with the infall mass profile. 



4. X-ray Data and Analysis 

A576 has been observed by both Einstein and ASCA. ASCA^s broad energy band 
(0.5-10.0 keV) is particularly useful for determining cluster temperatures. Because of the 
poor angular resolution of ASCA, we determine the emission weighted average temperature 
within 15', or ~ 0.5 /i~^Mpc. We obtained the screened data from GSFC. We extract a 
spectrum including all photons within a circle of radius 15' (60 pixels) centered on the 
cluster center for GIS data from a long (97 ksec) observation. The centroid of the SIS image 
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agrees with the position of the Einstein centroid within the ~ 0f4 uncertainty in the SIS 
position ( IGotthelf 1996| ). 



Using XSPEC (vlO.O), we fit the cluster spectrum to a model including absorption 
('wabs') parameterized by the column density of hydrogen (which we set to the galactic 
value) and the standard Rajmiond- Smith model ( [Raymond fc Smith 1977|) characterized 



by temperature, iron abundance, redshift, and a normalization factor. The iron abundance 
is measured relative to cosmic abundance. We fit the temperature, iron abundance, and 
normalization as free parameters. 

We use the weighting system developed by Churazov et al. (1996), which avoids 
rebinning the data into broad bins. Because there are few counts above 8.0 keV, we only 
include data from 0.8 8.0 keV, though we fit the spectrum to slightly different ranges to 
ensure that the fitted model parameters are consistent for different choices. 

We obtain acceptable fits by assuming that the gas is isothermal. More complicated 
models are thus unnecessary. Our best-fit model has an ICM temperature of 3.77±0.10 keV 
with an iron abundance 0.27±0.05 cosmic. Uncertainties are 68% confidence limits for one 



parameter. Figure |T0| shows the X-ray spectrum and the best-fit model. This temperature 
is 1.7(7 less than the temperature of 4.3±0.3 keV (David et al. 1993) from Einstein MFC 
data and 2.5a less than 4.02 ± 0.07 keV from an independent analysis of the ASCA data 
( [White 200q ). 

We calculate the expected flux between 0.01 and 100.0 keV from the best-fit model, 
yielding an essentially bolometric flux of 4.30 x 10~^^ ergs cm~^ s~^. We calculate a 
bolometric luminosity of 0.71 x lO^^/i"^ ergs s~^ (we assume go = 0.0) in agreement with 
0.73 X W^h-^ ergs s"^ from Einstem MFC data ( Pavid et al. 1993[) . The 0.3-3.5 keV 
luminosity is 0.46 x 10^'^/i~^ergs s~^, in excellent agreement with M96, who find a luminosity 
of 0.45 X 10^^/i~^ergs s~^ from the Einstein IFC data assuming an ICM temperature of 
4.3 keV. 

M96 found that a two-temperature Raymond-Smith model improved the fits to the 
Einstein spectrum. We test this result by adding a second thermal component to our 
isothermal model, but the fitting routine forces the temperature of the second component to 
be extremely small O.lkeV), suggesting that no second component is needed to explain 
the more complete ASCA data. This result contradicts the lower central temperature found 
with Einstein SSS data ( [Rothenflug et al. 1984[) . An independent anaysis of the ASCA data 



( [White 2000| ) suggests that A576 has a fiat temperature profile and no cooling flow. 

We use Einstein IFC data to analyze the surface brightness profile. Frevious studies 
( pones fc I* 'or man 1984] ; 1999; M96) fit the profile to the hydrostatic-isothermal j3x model 
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( lUavaliere fc Fusco-Femiano 1976 ) and found jS^ = 0.45 — 0.53, and 0.58-0.72 and core 
radius a = 50 — 70, and 90 — 160 /i~^kpc for Jones & Forman (1999) and M96 respectively. 
These models yield central gas and electron number densities of po = 3.6(2.3) x 10~^^ g 
cm~^ and 1.9(1.2) x 10~^cm~^ respectively. M96 suggest that the disagreement arises from 
fitting different radial ranges; Jones & Forman extend their analysis to larger radius than 
M96. We consider both models below. 



4.1. X-ray Mass Profile 

With the assumptions of hydrostatic equilibrium, negligible non-thermal pressure, and 
spherical symmetry, the gravitational mass inside a radius r is 

° jinipG ^ dhir dhir ' 



( [babricant, Lecar, fc Gorenstein 1980| ). For a uniform temperature distribution, the second 



term on the right hand side vanishes. We then only need to determine the gas temperature 
and the density distribution of the gas to calculate the gravitational mass. Under the 
standard hydrostatic-isothermal (3x model, the mass is related to (3x and a by 

3kT8 

where Mtot{< r) is the total gravitational mass within a radius r and the numerical 
approximation is valid for Tj^gy in keV and r and a in h~^Mpc. 

We use Equation ^ to calculate the total mass profile of A576 within 1.5h~^Mpc] the 
estimate outside ~ 0.5 h~^Mpc is an extrapolation. We use our 7456*^- derived temperature 
with the results of both M96 and Jones & Forman (1999) to estimate the mass profile. 
Figure |ll| shows these profiles as well as the gas mass profile for the M96 parameters. We 
also display an estimate of the cluster mass based on spherical deprojection (White, Jones 
& Forman 1997). This estimate is M ^ 1.07 x 10^^ H'^Mq at r = 0.299 /i^^Mpc and is 
larger than either (3^ model estimate. 

Evrard et al. (1996) develop a method to reduce the scatter in cluster mass estimates 
by relying solely on the emission-weighted gas temperature Tx within a radius where the 
mean density is 500 times the critical density. This radius, denoted by rsoo, varies with Tx 

as 

rsoo = (1.24±0.09)(^^)'^\-iMpc. (8) 
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The mass within r^QQ is approximately 

M5oo(Tx) = 1.11 X 10i5(^^)'^' h-'M^ (9) 

and has an average estimated-to-true mass ratio of 1.00 with a standard deviation of 8-15%. 
For A576, this procedure yields M500 = 2.57 x 10^^ H'^Mq at rgoo = 0.76 ± 0.16 h'^Mpc. 
This semi-empirical estimate of rsoo agrees with our more direct estimate of 
^500 = 0.96 ± 0.05 h^^Mpc from the caustics. From the infall mass profile, we 
estimate M500 = (5.1 ± 0.5) x 10^^ H'^Mq. 

Within ~ 0.5 h~^Mpc, the infall mass profile is a factor of ~2.5 larger than the 
isothermal X-ray mass profiles and a factor of ~1.8 larger than the deprojection estimate of 
White et al. (1997). M96 find a similar difference between X-ray mass estimates and other 
kinematic mass estimates of A576. These differences could be due to non-thermal pressure 
support, nonisothermality, or asymmetry. 



4.2. Gas Mass Fraction 

Using only X-ray data, the gas mass fraction fg (the fraction of total gravitational 



mass in hot gas) increases with radius (Figure |T2D, suggesting that the hot gas is more 



extended than the mass distribution, in agreement with previous studies ( Pavid, Jones, fc 



Forman 1995| ; [Markevitch fc Vikhlinin 1997| ; [Ettori fc Fabian 1999| ). The gas mass fraction 



is ~ 0.07/1 '^/^ at Rp = 0.5 h ^Mpc, less than (but consistent with) the average value of gas 
mass fractions in the most luminous clusters (e.g., [White et al. 1993| [Evrard 1997) ; [Mohr 



Mathiesen, fc Evrard 1999| ; [Ettori fc Fabian 1999| ). In hydrodynamical simulations (e.g. 



Evrard 1997] and references therein), the hot gas is more extended than the dark matter 



because of the shock heating originating during the merging of the halos which will form 
the cluster. In real systems, energy ejection from supernovae could also contribute to heat 
the gas. 

We observe a similar trend of increasing gas mass fraction with radius when we use the 
best-fit NFW infall mass profile to calculate the gas mass fraction (Figure |T2|). Because the 
infall mass profile exceeds the X-ray mass profile, the inferred gas mass fraction (~ 0.03/i^^/^ 
at Rp = 0.5 /i~^Mpc) is smaller. 
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5. Estimating the Contribution of Background Galaxies 

Within the completeness limit of our redshift survey, the elimination of background 
galaxies is straightforward; we simply remove non-member galaxies using the caustics in 
the projected radius-redshift diagram. We must estimate the background statistically at 
magnitudes fainter than mn = 16.5. For 16.5 < mn < 18.0, we only have spectroscopy in 
the central region (< 0.°75 radius; |Mohr et al. 2000]) . In this central region, the survey is 



91, 79, and 60% complete to limiting magnitudes of m/j = 17.0, 17.5, and 18.0 respectively. 
We assume that this region is sufficiently large to represent a fair sample of the background 
over the entire field. We then calculate the number of background galaxies in each 0.5 mag 
bin by assuming that the fraction of cluster members is the same for galaxies with and 
without redshifts. 

It is important to determine whether the survey of A576 contains any background 
or foreground groups or clusters of galaxies; these groups would lead to a localized 
enhancement in the number density of galaxies. Applying the adaptive kernel method 
( Pisani 1993| ) to all redshifts in our sample yields the estimated parent distribution shown 



in Figure |l^. A576 is very prominent. This distribution shows a background concentration 
of galaxies at 19, 700 km s~^; these galaxies concentrate in a region ~ 0.°50 — 0.°75 from 
the center of A576 (Figure P^). This concentration suggests an enhancement of background 
light in this region relative to a randomly selected field. Because we use this region to 
estimate the background, we may overestimate the background luminosity density across 
the entire region. However, the excess background may be present across the entire field. 

To estimate the variation of galaxy backgrounds in randomly selected fields, we analyze 
four 1° X 1° images of fields in the Century Survey (Brown et al. 200C) taken on the same 



nights with the same observing setup as the A576 data. Due to large-scale structure, the 
background counts vary more than expected from a Poisson distribution; the mean number 
of galaxies in a field is 121 and the variance is 19.2 (the four fields contain 149, 118, 110, 
and 107 galaxies). Weighted by luminosity, the variance is 21%. Figure ^ shows the 
magnitude distributions from these fields. For comparison, we show the estimate from the 
central region of A576. In the magnitude range of interest (16.5 < < 18.0), the estimate 
of background light from the average of the Century Survey fields is ~ 40% (roughly 2 a) 
smaller than that from the central region of A576. Because this enhanced background may 
or may not be present across the entire region, we estimate that the background subtraction 
is accurate to ~ 40% on scales ^ 1° square. 

A third method of estimating the background is to assume that the surface number 
density of cluster galaxies becomes negligible at the outermost radii of our survey; we 
require that the cluster number density be no less than that of known cluster galaxies. We 
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attribute additional surface number density to background galaxies. By measuring this 
number density as a function of limiting magnitude, we can estimate the number counts of 
background galaxies and the total background light. Figure ITBI shows the magnitude-number 



distribution from this method. Because this method assumes that no faint galaxies at large 
radii are infall members, this method should overestimate the background light. The ratio 
of the background light estimates calculated from the central region, the number density 
profile, and Century Survey fields are 1.14 : 1 : 0.67. 

The background estimates from the central region of A576 and from its outskirts agree 
surprisingly well. Thus, our determination of the background light seems to be a reasonable 
estimate of the background luminosity density across the entire region studied. We adopt 
the background values from the number density profile, which are intermediate between the 
values for the central region and the Century Survey. 

This estimate of the background is likely an overestimate. It is conservative in the sense 
that it tends to underestimate the cluster light. Because the amount of background light 
enclosed increases with radius, the effects of overestimating the background also increase 
with radius. 



6. Surface Number Density Profile 

One approach to studying the relative distributions of mass and light in a cluster is to 
compare the infall mass profile with the galaxy surface number density profile. We calculate 
the surface number density profiles with two different galaxy samples mig.s and mig where 
the subscript denotes the magnitude limit of the sample. For both samples, we exclude 
non- members with <16.5. The mie.s sample provides a lower limit on the asymptotic 
number density for fainter limits. We subtract a constant surface number density from the 
mis sample using the method described in §5. 

If light traces mass, the surface number density profile should closely resemble the infall 
mass profile (§3). The projection of the NFW profile yields a surface density profile S(-Rp) 

^j^Fw{R) = ^ — . [1 - X{R)] (10) 

^^"^^ ' 7rln(4/e)aW(i?2_i)^ ^ ^ ' 

where R = R/aNFW is the projected radius in units of the core radius, Nc is the number of 
galaxies within aNPWi and 

X{R) = (11) 
'R^ - 1 
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Similarly, the projection of the Hernquist profile can be written as 

27V 

^HerniR) = , , ^J {2 + R')X{R) - 3] (12) 

( IMahdavi et aTTOgp . 

Table |^ gives the results of our fits. For the NFW profile, the scale radius of the 
surface number density profile for the mie.s sample agrees with that for the mass profile. 
The scale radii for the mig sample and for the Hernquist fit to the mi6.5 sample are larger 
than the best-fit scale radii for the respective mass profiles. We display the surface number 
density profiles and the best-fit models in Figure There appears to be a relative excess 
of galaxies in the magnitude range 16.5 < niR < 18.0 between 1 and 2/;,~^Mpc. This excess 
leads to significantly larger scale radii for the mig sample than for the mig.s sample. This 
discrepancy may suggest the presence of a background cluster of galaxies (§5). The surface 
number density profile yields some evidence that galaxies are more extended than the dark 
matter, just as hot gas in the core is more extended than the dark matter. This evidence is 
stronger for the Hernquist model than for the NFW model; note that the Hernquist models 
produce better fits. One possible explanation of this difference is luminosity segregation, 
which we investigate next. 



7. Luminosity Segregation 

To convert to absolute magnitudes, we assume that all galaxies are at the distance of 
the cluster center; thus, Mr = uir — 35.29 — Ar + 5 log /i, where we estimate the dust 
extinction coefficient Aji from the dust maps of Schlegel, Finkbeiner, & Davis (1998) as 
listed in Table 1. We obtain the color excess at each galaxy position; the inferred extinction 
values range from 0.13 to 0.26 across the field. 



Figure shows the luminosity distribution for all galaxies with ran < 18 after removing 
all non-cluster galaxies with ttir <16.5, the completeness limit of our spectroscopic sample. 
The contribution of background galaxies steepens the distributions at Mr ^ —19. In the 
outer regions, the number counts begin to resemble N{'m) oc m 



0.6 



Luminosity segregation, if present, should be most apparent for the most luminous 
galaxies. In a large sample of clusters, Adami et al. (1998a) claim luminosity segregation 
among (on average) the four brightest members of each cluster. As a simple test of 
luminosity segregation. Figure |18| shows absolute magnitude versus radius for cluster 
galaxies in the complete mi6,5 sample. Significant segregation requires increasing absolute 
magnitude with radius. In fact, the three brightest galaxies in A576 are all > 1 /i~^Mpc 
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from the center, and the brightest cluster galaxy is 3.5 h~^Mpc from the center, more than 
two Abell radii away. The most luminous galaxy appears to be a normal SO galaxy, with 
an absorption dominated spectrum and no close companions. The next two most luminous 
galaxies also have absorption dominated spectra, but one of them has a close companion. 

As a quantitative test of luminosity segregation, we compare the cumulative luminosity 
distributions for all cluster galaxies with Mr < —19.0 (our spectroscopic completeness 
limit) inside and outside a radial cutoff (e.g.. Figure |1^). We vary this cutoff radius from 
0.05 — 2.50 h~^Mpc in steps of 0.05 h~^Mpc and perform a K-S test for each division. The 
only radial division that yields a separation at greater than a 90% confidence level is at 
1.45 h~^Mpc (95.2% confidence). However, the difference between the samples suggests 
luminosity anti-segregation, i.e., the outer sample is marginally brighter than the inner 
sample. 

To test the sensitivity of this result to the magnitude limit, we repeat the analysis 
using a magnitude cutoff of M/j < —20.7, approximately the value of M^, for the entire 
region (§8.1). With this cutoff, there is no significant difference at any radius (see Figure 
|19|). Combined with the result that the three most luminous galaxies in the sample are 
located far outside the core, we conclude that there is no significant evidence of luminosity 
segregation in A576. 



8. Luminosity Function and Light Profile 
8.1. Luminosity Function 



By correcting the luminosity function for the background counts (§5), we obtain the 



luminosity function shown in Figure |T^. Because we explicitly assume that few faint 
member galaxies are present at large radii, we are unable to constrain the faint-end slope of 
the luminosity function. 



We calculate the best-fit Schechter (1976) luminosity function. 



dn 



dM 



M)i 



(13) 



M 



in the range —22.5 < < —19, where is the characteristic absolute magnitude and a 
is the slope at the faint end. We find the best-fit form by minimizing x^- The parameters for 
the total luminosity function, M^, = — 20.7±0.4 and a = — 1.0±0.3 (95% confidence limits), 
are consistent with those of M96, LCRS ( [Lin et al. 1996| , a hybrid Gunn r-Kron-Cousins R 
system), and the Century Survey ( |Geller et al. 1997| , an R-selected survey). The fit yields 
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= 3.5 for 10 degrees of freedom; a single Schecter function is therefore a good model of 
the data. 



8.2. Light Profile 

We calculate the light profiles for the mig 5 and mig using the estimated background 
light from the asymptotic number counts around A576 (Figure ^0]). Using the absolute 
R-band magnitude of the Sun (M/? = 4.3 ; Zombeck 1990), we calculate the total light 
profiles for these subsamples. The potential overestimate of the background leads to an 
underestimate of the light profile, particularly at large radii. The mig sample yields an 
upper limit to the light profile by assuming that all galaxies without redshifts are cluster 
members (a lower limit comes from the rriR < 16.5 sample). 

Because the line-of-sight distribution of the galaxies is unknown, the observed light 
profile is the projection of the actual light profile. We therefore compare the observed light 
profile to the analytic forms in §5, replacing Nc with Lc, the amount of light enclosed within 
the scale radius a. The cumulative projected light profile for the Hernquist model has the 
simple form 

H< R,) = L,.,a^ ^/^^~/ (14) 
1 — it^ 

where Lfot is the total luminosity of the system ( Hernquist 1990| ). The results, shown in 



Table |^, reveal that the scale radii of the best-fit NFW profiles agree with that of the infall 
mass profile, whereas the Hernquist light profiles have scale radii a factor of ~ 1.8 — 3 larger 
than the best-fit mass profile. Figure ^ displays the best-fit profiles of each type. The 
values of rather large in all cases. It is evident from inspection of Figure I2T| that the 



observed light profile does not fall off as steeply as the NFW or Hernquist profiles at large 
radii {Rp ^ 2/i~^Mpc). If we convert the projected mass profile into a projected light profile 
by assuming a constant mass-to-light ratio (Figure ^OD, the predicted light profile increases 
more slowly than the observed light profile. This result suggests that the difference in scales 
between the mass and light profiles can not be explained by projection effects unless there 
are significant departures from spherical symmetry (e.g., our line of sight is perpendicular 
to the disk of an oblate system). 



9. Mass-to-Light Profile 

With the independently derived mass and light profiles, we can now calculate the 
R-band mass-to-light profile of A576 out to ~ 4 h'^Mpc, the limit of the infall mass profile. 
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By dividing the mass profile by the projected hght profile, we obtain the mass-to-light 
profile shown in Figure ^ M/Lr(< Rp) is largest in the core (~ 700 h) and decreases 
smoothly to the limiting radius of ~ 4 Mpc, where M/Lr ~ 300/i. We use the lower 
limits on the light profile to calculate "upper limits" on M/L/j(< Rp). These upper limits 
do not include the uncertainties in the mass profile; they demonstrate that the decreasing 
mass-to-light profile is not a result of underestimating the background light (in §5, we 
explain why we may overestimate the background light). We fit the profiles from the (mig) 
sample to a straight line (M/ Lr = aRp + b) and to a constant value (M/ Ln = b) using 
weighted least squares. Table ^ displays the results of these fits. The best-fit straight lines 
have negative slopes with high significance. Because the true M/L is always positive, the 
actual M/ L profile extrapolated to arbitrarily large radius cannot be a straight line with a 
negative slope. Again, because the values of M(< R) are not independent, the values of 
are only indicative. However, this exercise shows that a decreasing profile yields a much 
better fit than a fiat one. 

We note again that the light profile is the two-dimensional projection of the three- 
dimensional light profile; the deprojected mass-to-light profile of a spherically symmetric 



system would decrease more steeply than the ones presented here (see also Figure |20|) . As 
an example, we take the best-fit Hernquist mass profile, project it according to Equation 
(assuming spherical symmetry), and divide it by the projected light profile. The 



resulting projected- mass to projected-light profile, shown in Figure is larger in the inner 
~ l/i~^Mpc than the radial-mass to projected-light profile. 

This mass-to-light profile is an integrated profile, i.e., M/Lr{< Rp) = M(< Rp)/L{< 
Rp). The significance of the decreasing profile is more dramatic as a differential profile 
where M/L^^Rp) = 6M{Rp)/6L{Rp) (Figure |23|). Again, a decreasing mass-to-light profile 
yields a better fit than a flat one (Table |^). The outermost value of this proflle (~ 100 h) 
should be closest to an estimate of the universal value of M/Lr. Interestingly, this value 
agrees with the mass-to-light ratios of some elliptical galaxies ( |Mushotzky et al. 1994| ; 
Bahcall, Lubin, fc Dorman 1995). 



10. Discussion 

A decreasing mass-to-light proflle is perhaps surprising. One immediately wonders what 
systematic effects could mimic a decreasing M/Lr proflle. Is the mass proflle incorrect? 
The infall mass proflle yields larger values than X-ray estimates in the inner regions and 
lower values than virial theorem estimates in the outer region. Non-thermal pressure 
support, nonisothermality, and deviations from virial equilibrium are plausible explanations 
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for these effects. Another unknown systematic effect is the actual distribution of galaxy 
orbits. However, simulations suggest that this effect is not large (D99). The infall mass 
profile itself is uncertain by a factor of 2. Hidden systematic effects could affect the shape of 
the mass profile, a possibility we plan to explore with other systems. The infall mass profile 
is unstable from 1.5 — 2.2h~^ Mpc, but beyond this range, M/Lr continues to decrease 
even though the mass profile is well determined. A Jeans analysis of the virial region will 
provide an interesting check on the infall mass profile ( pVIohr et al. 200d| ); preliminary results 
suggest good agreement between the infall mass profile and the Jeans mass profile. We note 
also that the X-ray mass found by White et al. (1997) provides a lower mass limit at small 
radius. While this mass estimate is smaller than the infall mass estimate, it is not sufficient 
to reconcile the observations with a constant mass-to-light ratio. 

Is the light profile poorly estimated? In the inner regions, we use SExtractor 
magnitudes, which may underestimate the diffuse light associated with Brightest Cluster 
Galaxies (BCGs) by ~ 50% ([Gonzalez et al. 2000|) . Even though A576 has no cD galaxy. 



diffuse light near the center of the cluster may partially account for the observed decreasing 
M/L profile. The magnitude of this effect (~ 20% of the total light within 1.5 /i~^Mpc), 
however, is not sufficient to reconcile the observed M/L^ profile with a constant value. 
Background subtraction always contains some uncertainty; we show the limits of this 
uncertainty based on galaxies reliably contained within the infall region. The background 
light we assume exceeds that in randomly selected fields, and two independent estimates of 
the background agree remarkably well. The shape of the MjL^ profile is insensitive to the 
magnitude limit. The surface number density profile also supports a larger scale for the 
galaxies than for the mass. We note that any systematic underestimate of total galaxy light 
(e.g., missing halo light) would likely affect the magnitudes of all galaxies. Such an effect 
would change the absolute value of M/L/j but would be unlikely to alter the shape of the 
profile. 

Are projection effects skewing our results? According to the infall model, the mass 
profile is a true radial profile; the light profile is the two-dimensional projection of the 
actual light profile. Assuming that luminosity density decreases monotonically with radius, 
the projected light profile is steeper than the actual light profile. This effect would increase 
the significance of the decreasing mass-to-light profiles presented here. Departures from 
spherical symmetry could also affect the mass-to-light profile; looking down the barrel of a 
prolate spheroid would exaggerate the decrease in M/L; looking at the disc of an oblate 
spheroid would mitigate it. There are no obvious indications of oblate or prolate structure 
in the sky distribution of infall members (Figure |^). Analysis of other systems should 
provide an important check on these projection effects. 
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Is the decreasing mass-to-light profile an effect of calculating the luminosity from 
R-band photometry? Galaxies in the cores of clusters are redder than those far outside 
the cores. This effect causes the R-band light density profile to decrease more steeply with 
radius than a B-band light density profile. This effect would oppose the apparent trend of 
a radially decreasing mass-to-light profile. 

M96 demonstrate that the emission-dominated galaxies have different kinematic 
properties than absorption-dominated galaxies. On this basis, M96 removes the light 
contribution from emission-dominated galaxies in calculating the mass-to-light ratio. 
Because many of these emission-dominated galaxies are in the infall region but probably 
not inside the virial radius, the mass-to-light profile is biased low in the central regions and 
high in the outer regions. 

Finally, we note a curious feature of the D99 simulations. For ACDM models, the 
mass-to-light profiles (in B band relative to the global value) show a decreasing trend similar 
to that shown here. Clusters in a tCDM model produce the opposite effect, mass-to-light 
profiles which increase with radius. D99 attributes this effect to the deficit of blue galaxies 
in the centers of clusters in the ACDM model, but later analysis reveals a similar effect 
for mass-to-light profiles measured in I band. This result suggests that the decreasing 
mass-to-light profile may contain information on the underlying cosmology. 



11. Conclusions 



We calculate the mass profile of Abell 576 to ~ 4 h^^Mpc from the observed infall 
pattern in redshift space. The amplitude of the resulting mass profile is larger than 
determined from X-ray observations to their limiting radius of 0.5 Mpc and smaller 
than virial mass estimates at larger radius. The infall mass profile agrees extremely well 
with an NFW profile or a Hernquist (1990) profile, and it is strongly inconsistent with an 
isothermal sphere profile. This result agrees well with a similar analysis of Coma by GDK. 
Our best-fit mass profile implies that the fraction of gravitational mass contained in hot 
gas increases with radius to the limit of the X-ray data, in agreement with previous studies 
of other clusters (e.g., [David, Jones, fc Ibrman 1995|; [Markevitch fc Vikhlinin 1997]; [Ettori 



fc I'abian 1999[ ) and simulations ( | White et al. 199B| ; [Evrard 1997] ). The hot gas therefore 
appears to be more extended than the dark matter. 

The decreasing amplitude of the caustics suggests that the mass increases more slowly 
than for an isothermal sphere. This inference is supported by fitting mass profiles to 
the data. GDK found that the NFW profile describes the mass profile of Coma much 
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more accurately than an isothermal sphere profile, and Figure 1 of CYE clearly shows the 
existence of radially decreasing caustics in the composite CNOC cluster. Many clusters 
therefore show evidence of mass profiles shallower than an isothermal sphere; the masses of 
clusters probably increase no faster than \n{Rp) at large radius. 

Using photometric data from a large area surrounding the cluster, we find little evidence 
for luminosity segregation. In fact, there is some evidence of luminosity antisegregation; the 
three brightest cluster galaxies lie far outside the center of the cluster. 

The R-band mass-to-light ratio is largest in the core (~ 700 h) and decreases 
smoothly to the limiting radius ~ 4 Mpc, where M/Lji{< Rp) ~ 300/i. The 
decrease is more dramatic in the differential mass-to-light profile. At the limit of our 
sample, dM/dL ~ 50 — ISO/i. The luminosity density from the Century Survey is 
PL = 0*L*r(2 + a) = 6.0 ± 2.6 x lO^h^ Mpc'^ which yields {M/L)crit = 459 ± 199/i ([GellCT 
et al. 199^ ). Our estimate of the global value of M/L thus implies Qm ^ 0.22 ± 0.1, or 

^ 0.4 at a 95% confidence level. This estimate of flm is an upper hmit due to the 
unknown contribution of faint {rriji > 18) galaxies. Note that using the luminosity function 
parameters from the LCRS ( |Lin et al. 1996| ) yields a significantly smaller luminosity density 
corresponding to {M/L)crit = 1000 ± 150h. Adopting this luminosity density reduces our 
estimate of the matter density to Qm ^ 0.10 ± 0.05. The decreasing M/L profile suggests 
that the dark matter is more concentrated than the optical light of the galaxies. We 
consider possible systematic effects (§10) and find that they generally oppose the decreasing 
mass-to-light ratio. If general, our results place strong constraints on possible variations of 
mass-to-light ratios with scale as well as on the global value of Qm- 

Photometric observations in other bandpasses would allow an examination of the 
effects of stellar populations on the mass-to-light profile. Similar studies of other clusters, 
particularly studies of clusters with better defined infall regions at large radius, should test 
the generality of our results. 

This project would not have been possible without the assistance of Perry Berlind 
and Michael Calkins, the remote observers at FLWO, and Susan Tokarz, who processed 
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the Chandra Science Center. The Chandra Science Center is operated by the Smithsonian 
Astrophysical Observatory for NASA under contract NAS8-39073. AD was supported by a 
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Fig. 1. — Comparison of magnitudes from the two nights for galaxies in the central region. 
For rriji < 18, the zero-point offset is -0.002 mag with a ^ 0.09 mag. The scatter increases 
with apparent magnitude. 
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Fig. 2. — (a) Comparison of our MAG_BEST magnitudes with those of M96. Our magnitudes 
are systematically brighter, (b) Comparison of our magnitudes with those of M96 after 
processing both data sets with SExtractor. 
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Fig. 3. — Fraction of galaxies with measured redshifts as a function of mu magnitude. The 
sohd hne shows galaxies with Rp < 1°, the dashed hne shows galaxies at larger radii. We 
also display the total number of galaxies in each bin (thick lines). 
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Fig. 4. — Redshifts as a function of projected radius in A576. The solid lines are the caustics 
determined from our adaptive kernel estimate with 1-a error bars. 
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Fig. 5. — Dependence of caustic location on the parameter g, the ratio of velocity uncertainty 
to positional uncertainty. From top to bottom, the caustics are fit with q = 10, 25, and 50. 
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Fig. 6. — Distribution of cluster and field galaxies on the sky. Open circles, filled circles, 
and crosses represent cluster members, members with rriR < 16.5, and background galaxies 
respectively. 
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A576 




Rp (h-i Mpc) 

Fig. 7. — Mass profile of A576. Tlie filled squares are from the adaptive kernel estimate of 
the caustics with 1 — a uncertainties shown. The upper and lower filled hexagons are virial 
mass estimates (Girardi et al. 1998) omitting and including the surface term respectively; 
the filled hexagon at small radius is an interpolation by Girardi et al. discussed in the text. 
Crosses show our virial and projected mass estimates, the star shows an X-ray mass estimate 
(§4.1, [White, Jones fc Forman 1997|) . The solid, dashed, and dash-dot lines are the best-fit 
NFW, Hernquist, and isothermal sphere profiles respectively. 
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Fig. 8. — Mass density profile of A576. Lines have the same definitions as in Figure |^ 
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Fig. 9. — Velocity dispersion profile for A576. Open squares show the velocity dispersion 
profile ap{Rp). Filled squares show the cumulative velocity dispersion profile (Tp{< Rp). 
Dashed and dash-dotted lines indicate the predicted profiles for the best-fit Hernquist and 
NFW profiles respectively assuming isotropic orbits (/? = 0). 




Fig. 10. — ASCA GIS spectrum of A576. The solid line is the best-fit single-temperature 
Raymond-Smith thermal plasma model. 
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Fig. 11. — Mass profile of A576 derived from X-ray observations under the assumption of the 
hydrostatic-isothermal model and parameters from Jones & Forman (1999, dash-dotted 
hne) and M96 (dashed line). The sohd hne is the gas mass profile and the vertical line at 
Rp — 0.5 /i~^Mpc shows the hmit of the X-ray data. Filled squares show the infall mass 
profile. The open hexagon and star are X-ray mass estimates from White et al. (1997) and 
from the estimator of Evrard et al. (1996) respectively. 




Fig. 12. — Gas mass fraction profile. Mgas is taken from the X-ray temperature and density 
profile, Mtot is the best-fit X-ray mass profile (dashed line) or Hernquist (solid) or NFW 
(dash-dot) infall mass profile. The vertical line is the limit of the X-ray data. 



-42 - 




Fig. 13. — Adaptive kernel estimate of fka{cz), the velocity distribution function. Arrows 
indicate the peak and limits of A576 and the limits of a background concentration of galaxies. 
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Fig. 14. — Distribution on the sky of background concentration of galaxies. The cross marks 
the X-ray center of A576. 



14 15 16 17 18 19 

Fig. 15. — Apparent magnitude distribution of four randomly selected fields from the Century 
Survey (solid lines). The dashed line is the background estimated from the central region of 
A576; the dash-dot line is estimated from the asymptotic number density profile. 
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Fig. 16. — Surface number density profile for A576. Squares show the profile calculated from 
niR < 16.5, hexagons show the background-subtracted profile for uir < 18. The sohd and 
dashed lines are the best fit NFW and Hernquist profiles respectively for the niR < 16.5 
sample. 



Luminosity Function 




Fig. 17. — Luminosity function for entire region. Filled and open squares are calculated 
with and without background subtraction respectively. The vertical line is the limit of our 
survey. The curve indicates the best-fit Schechter luminosity function for Mr < —19. 
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Fig. 18. — Absolute magnitude versus projected radius. Significant luminosity segregation 
requires that absolute magnitudes increase with radius. The three most luminous galaxies 
all have Rp > 1 /i^^Mpc. 
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Fig. 19. — (a) Cumulative absolute magnitude distributions of galaxies inside (solid lines) 
and outside (dashed lines) Rdiv = 1-45 h^^Mpc. The outer sample is brighter than the inner 
with 95% confidence for Mr < -19.0. (b) Same for Raiv = 1.45 h-^Mpc and Mr < -20.7. 



Table 1. Photometric and Spectroscopic Data^ 



RA DEC cz acz tur Ar 

(J2000) (J2000) (kms-^) (km s"^) 



7 10 58.32 


+56 29 25.4 


39154 


37 


15.84 


0.06 


0.15 


7 10 59.40 


+56 30 15.8 


14569 


17 


16.11 


0.06 


0.15 


7 11 03.94 


+56 39 38.2 


16042 


15 


16.38 


0.06 


0.13 


7 11 04.13 


+56 39 46.4 


15918 


16 


15.92 


0.06 


0.13 


7 11 09.26 


+56 02 13.2 


13780 


29 


15.59 


0.06 


0.16 



^The complete version of this table is in the electronic edition of the 
Journal. The printed edition contains only a sample. 
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Fig. 20. — R-band luminosity profile of A576 (crosses). We exclude all non-cluster galaxies 
with rriR < 16.5. The faint solid and dash-dot lines are the profiles calculated from rriR < 18 
galaxies with and without background subtraction respectively. The heavy solid and dashed 
lines are the best-fit Hernquist light and mass profiles (converted with an arbitrary mass-to- 
hght ratio). 
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Fig. 21. — R-band luminosity density profile of A576. Symbols are as in Figure 16. The solid 
and dashed lines are the best fit NFW and Hernquist profiles respectively for the rriR < 16.5 
sample. 
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Fig. 22. — Cumulative mass-to-light profile of A576 calculated from the infall mass profile 
and the projected light profile. Squares estimate the light profile from galaxies with thr < 18; 
upper limits include only light from confirmed members {rriR < 16.5). The dash-dotted line 
indicates the projected best-fit Hernquist mass profile divided by the projected light profile. 
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Fig. 23. — Differential mass-to-liglit profile. Squares represent the mis luminosity profile. 
Upper limits reflect the minimum light associated with the cluster. 
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Table 2. Spectroscopic Data^ 



RA 


DEC 


cz 




(J2000) 


(J2000) 


(km s-^) 


(km s"''') 



6 53 01.56 


+56 15 35.3 


10057 


32 


6 54 33.82 


+55 14 41.6 


14215 


39 


6 55 15.14 


+55 28 42.6 


16899 


29 


6 55 48.17 


+56 04 56.3 


19347 


43 


6 56 00.38 


+55 24 02.2 


7908 


26 



'^The complete version of tliis table is in the 
electronic edition of the Journal. The printed 
edition contains only a sample. 



Table 3: Mass Profile Fit Parameters 



Profile 


a 


95% 


C 


95% 




V 




/i~^Mpc 


/i~^Mpc 


lO^^Mo/Mpc 


lO^^Mo/Mpc 






NFW 


0.13 


0.11-0.16 


2.9 


2.7-3.2 


7.23 


25 


Hernquist 


0.43 


0.37-0.49 


2.5 


2.3-2.7 


1.64 


25 


Isothermal 






0.322 


0.308-0.334 


179 


26 



Table 4: Surfa ce Number Density Profile 

Profile Sample a 95% ^ Prob 

NFW mi6.5 124 ol2 0.08-0.18 33J 20 (103" 

NFW mi8 340 0.42 0.24-0.74 10.4 20 0.96 

Hernquist wie.s 605 0.68 0.58-0.82 28.9 20 0.09 

Hernquist mis 1060 1.20 0.80-1.64 11.5 20 0.93 
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Table 5: Projected Light Profile 



Profile 


Sample 


a 


95% 


Lc 


95% 




V 






/i~^Mpc 


/i~^Mpc 


IQiOio 


lO^OLo 






NFW 


mi6.5 


0.15 


0.12-0.18 


59 


56-62 


333 


22 


NFW 


mi8 


0.45 


0.34-0.59 


152 


124-188 


17.9 


22 


Hernquist 


TO16.5 


0.76 


0.70-0.84 


260 


252-272 


347 


22 


Hernquist 


mi8 


1.27 


1.03-1.59 


470 


400-560 


26.8 


22 


CumHern 


mi6.5 


0.86 


0.85-0.87 


245 


243-248 


2524 


365 


CumHern 


mi8 


1.47 


1.45-1.48 


459 


456-463 


6121 


1846 



Table 6: Fits of M/L Profile to MjLn = aRp + b 



Type 


Sample 


a 


b 


x' 


1/ 


Int 


mie.5 


-62±4 


651±10 


2.0 


22 


Int 


mis 


-90±5 


592±13 


5.5 


22 


Int 


mi6.5 




516±14 


25 


23 


Int 


mi8 




378±10 


81 


23 


Diff 


mi6.5 


-158±31 


583±67 


127 


22 


DifT 


Wis 


-109±27 


376±61 


206 


22 


Diff 


TO16.5 




263±10 


278 


23 


Diff 


mi8 




145±6 


354 


23 



